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CHAPTER 1 


INTRODUCTION 


Beam-like structures cantilevered to the Shuttle cargo bay will, 
in general, experience axial loads as a result o£ attitude maneuvers. 

In our previous \»ork (Reference 1), we developed dynamics models for 
such a class of orbiter-payload systems with the effects of axial loads 
ignored for reasons of simplicity and tractability. In this report, we 
examine that approximation in a more simple setting which clearly ex- 
poses the phenomena of interest. Presented are analyses of the planar 
transverse bendirg behavior of a uniform cantilevered beam with rigid 
tip body under the action of constant axial base acceleration. Exact 
steady-state and free vibration solutions are obtained for various 
forms of this problem with additional approximate solutions for selec- 
ted cases. While the results of this study are of limited direct ap- 
plicability to the problem of Reference 1, they are of general inter- 
est, provide much useful insight, and serve as a basis from which to 
address more complicated situations. 

The beam i •> taken to be long, slender, and inextensible with 
uniform mass, stiffness, and cross section. Only small elastic trans- 
verse bending deformations in a plane of symmetry of the cross section 
are considered. The tip body is rigidly attached to the beam, with its 
mass center having an arbitrary offset with respect to the attachment 
point. In all cases, the base acceleration is constant in magnitude 
and directed along the undeformed longitudinal axis of the beam. The 
governing partial differential equation is shown to be linear with 
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variable coefficients, accompanied by nonhomogeneous boundary condi- 
tions for the general problem of arbitrary mass center offset. The 
steady-state response is examined for a number of beam end conditions. 
For the cases of free-end, tip mass, and tip body with mass center 
along the beam tip tangent line (henceforth referred to as restricted 
mass center offset), the boundary value problems are homogeneous. 

These classical Euler buckling problems are solved e cactly for the 
critical buckling loads/accelerations. For the problem of arbitrary 
tip body mass center offset, it is shown that a unique steady-state 
solution exists, except for certain critical values of the base accel- 
eration for which no such solution exists. The treatment of transverse 
vibration begins with the case of restricted mass ce. ter offset. Tlie 
boundary value problem is homogeneous and admits the usual separation 
of variables solution with harmonic time dependence. Uie ensuing vari- 
able coefficient ordinary differential equation for the spatial func- 
tions is solved exactly in terms of a power series. Application of the 
boundary conditions yields the characteristic equaticn and eigenfunc- 
tions. For the case of arbitrary tip body mass center offset, the 
vibration problem is nonhomogeneous. It is shown thai the exact solu- 
tion may be written as the sum of the steady-state solution obtained 
previously, and a superposition of simple harmonic mccions which cor- 
respond to the exact solution of the associated homogeneous problem. 
This approach, which takes advantage of the steady-state solution, 
yields a greatly simplified form of the final solution. An assumed 
modes formulation is detailed for the restricted mass center offset 
case. Ihe approximate solutions to the free vibratioi. and buckling 
characteristics serve to check the exact analyses. A particularly 
useful and general method for recovering the strain energy from the 
governing partial differential equation and boundary conditions is set 
forth. Selected FORTRAN programs and numerical results are provided. 
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CHAPTER 2 


MOTION EQUATIONS 


2.1 Derivation of the Partial Differential Equation 
and Boundary Conditions 


In this chapter, we derive the partial differential equation and 
boundary conditions governing the planar transverse vibration of a beam 
with tip body under the action of an axial force. Figure 2-1 depicts 
the situation in the undeformed configuration. At one end of the beam, 
(x = 0) , an axial force is applied imparting a velocity, vgCt), along 
the x*axis. The beam is assumed inextensibls and of length 1. At the 
other end of the beam, (x = £), a tip body is rigidly attached at point 
P. Hie rigid tip body has mass m and moment of inertia about its mass 
center I^-. Hie distance from P to the mass center of the tip IxxJy is 
c, and this directed line segment makes an angle, Y, with the beam tip 
tangent at P. Figure 2-2 shows the system in the deformed state. The 
base 0 of the beam moves along the x direction while the deformation 
u(x, t) , occurs along the y direction. Note that the angle Y is main- 
tained constant. 




I 

INERTIAL X 

FRAME 



Figure 2-1 . 
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Figure 2-2. 

Figure 2-3 is a free-body diagram of a beam element between x 
and X + Ax. T(x, t) is the internal tension, s(x, t) is the shear 
force, and M(x, t) is the bending moment. Since the beam is assumed to 
have no deformation along the x axis, the acceleration component along 
X IS aQ(t) for any point in the beam, where ^^(t) is the prescribed 
base acceleration at 0. 


s(x + ax, t) 



Figure 2-3. 
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Equilibrium of the beam element along the x direction gives 


3t 


( 2 - 1 ) 


where p is the mass per unit length of the beam. 

Similarly, equilibrium along the y-direction gives 


9S 3^u 


( 2 - 2 ) 


Applying rotational equilibrium to the beam element and neglecting 
rotary inertia we obtain 


3M 3u 

^ + S(x,t) - T(x,t) 37 = 0 


(2-3) 


Differentiating Eq. (2-3) and using Eq. (2-2) , we obtain 


3^M 3^u 3 r 3ui - 

3x 3t 


Assuming an Euler-Bernoulli beam 


M(x, t) = EI(x) 


Iji 

3x2 


we arrive at the partial differential equation for the beam deflection 


3x 3x 3t 


(2-4) 
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c* 


We also have the differential relations for the natural boundary condi- 
tions 


M(x,t) = EI(x) 


afu 

3x^ 


(2-5) 


S(x,t) = 


2 

9u 9 r 9 un 

T(x,t) [ei(x) — -] 


9x 


( 2 - 6 ) 


Kinematics of the Tip Body 


Let Vp and denote the inertial velocity of P and the angular 
velocity of the tip body, respectively 


V 

P 


v^(t)t 



(^ft) 3 


( 2 - 7 ) 


-*• 

0) 

p 


95i_ 

9t9x 


(Jl,t) k 


( 2 - 8 ) 


^ 

where i, j, k denote unit vectors along the directions of x, y, z, re- 
spectively. 

If we let c be the vector from P to the mass center of the tip body and 
-► 

v^ the velocity of the tip body mass center then 


= v_ + 0 )_ X c 
P P 




- c sin Y 


9 u 
9t9x 


(Z,t)]i + [|^ [l,t) + c cos Y 


^2 
9 u 

9t9x 


(^f t) ] ] 

( 2 - 9 ) 
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+ <0 


3u 3^u 

where we have dropped the nonlinear term (Jl,t) -r— r- (£,t) . 

oX Oudx 

Differentiating Eqs. (2-8) and (2-9 ) , we have the acceleration expres- 
sions 


3 u 
3t^3x 


a.t) it 


( 2 - 10 ) 


[ao(t) 


o3 
3 u 


.3 2 

c sin T — :^(Jl,t)]i + [ — ^(A,t) + c cos y 

3t 3t 3x 


3t^3x 


(^,t)] 

( 2 - 11 ) 


Natural Boundary Conditions at x = £ 


Figure 2-4 is a free body diagram of the tip body indicating the 
force and moment exerted by the beam at P. 



Figure 2-4. 
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eft 


The equation of motion of the tip body when resolved along the 
directions of x and y yields the two scalar equations 


9 u 

-T(Jl,t) = ma (t) - roc sin Y 

° 9x9t^ 


( 2 - 12 ) 


g2 g3 

-S(£,t) = m — ^ (A,t) + roc cos Y ^ (^rt) 


(2-13) 


at 


9x9t 


Using Eqs. (2-6) and (2-12) allows us to write the first boundary 
condition at x = ?. as 


g3 .3 

roc cos Y ^ (^•/t) - El — ^ (^.,t) 


9x9t 


3x" 




g2 

+ m — ^ (£,t) 
3t 


= 0 


(2-14) 


where we have assumed a uniform El and have dropped the nonlinear term 
3 3 ^ 

me sin Y (£ft) ^ (£»t). 

3x3t^ 


It can be shown that the z component of the time rate of change 
of the angular monentum of the tip body about its mass center is given 
by 


- 

^ 3x3t^ 


(Jt,t) 


If is the net moment about the mass center of the tip body 

5 = -M(£,t)k + c X s(£,t)3 + c X T(£,t)i 

© 
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Using Eqs. (2-5) and (2-6) and dropping nonlinear terms we obtain 




g2 3 ^ 

= -[ei — ^ (Jlft) + El c cos Y — ^ {H,t) + c sin y T(il,t) ]k 


3x 


3x 


Using Eq. (2-12) 
at X = £ is 

2 2 

(I^ + me sin 


(2-15) 

for T(£,t) in Eq. (2-15), the second boundary condition 

g3 3 2 

Y) ^ (A,t) + EI c cos Y — ^ (£,t) + EI (£,t) 

3x3t 3x 3x 

- me sin y a {t) =0 (2-16) 

o 


The partial differential equation governing u(x,t), Eq. (2-4), 
requires a specification of the internal tension T(x,t) for 0 ^ x ^ 
Integrating Eq. (2-1) and using Eq. (2-12) for a boundary condition we 
obtain 


33 

T(x,t) = -[p(£ - x) + m]a (t) + me sin y ^ (£,t) 

° 3x3t^ 


(2-17) 


Assume we have a compressive load, PQ(t) > 0 applied at x = 0. Then 


T(0, t) = - P^(t) 


Hence 


P (t) 
o 


(p£ + m)a (t) - me sin y 
o 


3^u 

3x3t^ 


(£,t) 


(2-18) 


which relates the applied axial force to the base acceleration. 
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The governing equation of motion, Eq. (2-4) , can be rewritten 
with the aid of Eqs. (2-17) and (2-18) as 


El 


° 


(t) [(1 - 


px 


p£ + m 


9u -| 


^2 
O u 

P — = 0 

at 


(2-19) 


We assume the beam to be cantilevered at x = 0 so that 

3u 

u(0,t) = -^ (0, t) =0 for all t (2-20) 

The partial differential Eq. (2-19) ’ s to be solved subject to 
the geometric boundary conditions of Eq. (2 -20) and the natural boundary 
conditions of Eqs. (2-14) and (2-16). 

In Reference 1, we studied the transverse vibration of a clamped 
beam with tip body with no base acceleration. On page 10 of that report 
it was stated that y must be restricted to ?ero in order that no axial 
loads be introduced into the problem. With the present analysis, one 
can accommodate any axial loads that may be introduced by y ^0, and 
will discover, surprisingly, that the results given in Reference 1 can 
be directly extended to the case of y * 0. 
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CHAPTER 3 


STEADY-S t’ATE SOLUTIONS 


3.1 Buckling of a Cantilevered-Free Beam Under Axial Thrust: 

In this section, the possibility of a cantilevered-free beam (no 
tip body) buckling under axial thiust is investigated. The governing 
eguation is Eg. (2-19) with no time dependence and m = 0. 


El 


dx'^ 


+ P 


cr 


I; f’ - SI - 


dy-i 


dx ' 


(3-1) 


where P denotes those values of thrust which lead to buckling, 
cr 

The two geometric boundary conditions at x = 0 are 


y(0) = y'(0) = 0 

The boundary conditions at x = £ are obtained from Eqs. (2-14) and 
(2-16) with m = c = 0 


y'*(£) = y'"(£) = 0 

Integrating Eg. (3-1) and setting the constant of integration to zero we 
obtain 


^ <1 - kA) S . 0 

, 3 El dx 

dx 


(3-2) 
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subject to the three boundary conditions 


Y(0) = y'(0) = y"{£) = 0 


making the substitution 


^ - 


dx 


w, 


■ M'" 


(1 - x/£) 


transforms Eq. ( 3 - 2 ) into 


d^w 

dz^ 


+ zw = 0 


with boundary conditions 


w = 0 at z = 


'p a 
cr 

^ El 


2 \ 1/3 


— = 0 at z = 0 

dz 


The general solution of the differential equation can be written 


where c^ and C2 are arbitrary constants and J1/3 denotes the Bes- 
sel function of the first kind of order 1 / 3 . 
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0 . 


The boundary condition — = 0 at z = 0 implies that c, = 

dz 1 

Let denote the roots of the Bessel function of the first kind of 
order (-1/3) » i.e., 


^-1/3 '^n> 


0 


Applying the remaining boundary condition we find that the critical 
buckling loads are given by 


P 

cr ,n 


9 .2 El 
4 


n — 1f2r3f 


(3-3) 


The first three buckling loads are 


cri 


El 

7.8664 


cr2 


El 

55.977 

Z^ 


cr3 


El 

148.51 


3.2 Buckling of Cantilevered Beam with Tip Mass 
Onder Axial Thrust 

The governing equation is Eq. (2-19) with no time dependence and 
the boundary conditions ate Eqs. (2-20), (2-14), and (2-16), with 1^=0 
and c = 0. 


d^y 

El — f + 
dx 


"cr dx 


[(1 - 


px 


-1 = 
pZ + m-^ dx-* 


(3-4) 


mP 

y(0) = y' (0) = 0 : y”(Z) = o ; ei y"'(Ji) + y' (^) = 0 

where we have used Eq. (2-18) to express Sq in terms of Per* 
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Integrating Eq. (3-4) and setting the constant of integration to 
zero, we will satisfy the last boundary condition. 


dx 


Z + iZ = 0 

3 El ^ pi. + 10 -^ dx 


If we make the substitutions = w and 

dx 


dy 

= w 


r cr cpJi +m'|2i1/3/', px 
“ 1-Ei 0 J i “ p£ + m-* 


"El ^ p 


this equation transforms to 


d^w 

dz2 


+ zw = 0 


(3-5) 


with boundary conditions 


w = 0 at z = Zi 


and 


dw 

dz 


0 at z 


z 


2 


where 


/'^cr-,1/3 r P^ + m -> 2/3 
^EI ^ ^ P •' 

and 

, 5L r__!£L_ll/3 

p^^-^ EI(p£ + m) 
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Pg has already been shown 


w 


dw 

dz 




/3V3 “ ^ °2 


z [c J f— Z^^^} - c 

' 1 -2/3 3 2 




J (- z 

2/3^3 


3/2 „ 


Applying the boundary condxtions, we obtain a homogeneous system in 
and c^. Setting the determinant to zero (we want y(x) ^0), we obtain 
the following transcendental equation for the critical buckling loads. 


J jfd . J (. 


m 




/1 + 1/m* 

+ * -*U) • J_2,3 (■ 


/l + 1/m* 


Xl = 0 (3-6) 


where 


m* = m/p£ 


and 


X 



Note that the roots of fK[. (3-6) depend only upon the ratio of the tip 
ndss to the beam mass. 

In the limit as m* -*-0, the critical buckling loads in Eq. (3-6) 
should approach those found previously for a cantilevered- free beam 
(Rq. 3-3). This is readily established if we recall the facts 

lim J_,,(x) = 0 

X+0 2/3 

and 

lim J ^ (x) = 00 

x^ -2/3 
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3.3 Buckling of Cantilevered Beam with Tip Body 
Under Axial Thrust (y = 0) 


The governing equation is Eq. (2-19) with no time dependence and 
the boundary conditions are Egs. (2-14), (2-16), and (2-20) with Y = 0. 


El . p I- [ ' - *1 = 0 

,4 cr dx ' p£ + m' dx-' 
dx 


y(0) = y' (0) = 0 ; El c y' ’ ' (£) + El y" (£) = 0 (3-6) 


El y" • (£) + n. 


cr 


p£ + m 


y'(£) = 0 


(3-7) 


Integrating the differential equation and setting the constant of 
integration to zero, we will satisfy the boundary condition Eq. (3-7). 


P - 

d_y _c£ c px >, ^ 

^3 El '■ p£ + dx 
dx 


= 0 


If we make the same substitutions as in the previous section we obtain 


d^w 

dz^ 


+ zw = 0; w = 0 at z = z. 


Substituting Eq. (3-7) into Eq. (3-0) allows us to write the second 
boundary condition as 

(p£ + m) El y"(£) - me P y' (Ji) = 0 (3-8) 

cr 


or in terms of w(z) 

[(p£ + m)EI]^/^ (P^^p)^/^ (z_) + me P^^w(z_) = 0 (3-9) 

cr dz 2 cr 2 
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Using the results from the previous section. Eg. (3-9) can be written as 


>^1 + ra* J ( 


m’ 


2/3 . , : X) + I C* X J ( X) 

/I + 1/m* ^ /I + 1/m* 


li /i? c* X J , ,, f 2^ 

r ■ /I + 1/m* 

I 


- •'1 + m* J2/3 ^ 

/I + 1/m* 


= 2 = ° 


(3-10) 


where c* = c/£. 

The first boundary condition at z = zi can be written as 

J,/3[(1 + m*)X)ci + J-i/3((1 + m*)X)c2 = 0 (3-11) 

The homogeneous system, Eqs. (3-10) and (3-11), will have a nontrivial 
solution if 


J-1/3(aX) [>^ J_2^3(3X) + -| c* X J^^^(eX)] 

+ Jl/ 3 (aX) [✓a J 2 ^ 3 ( 0 X) - I c* X J-i/ 3 (eX) ] = 0 (3-12) 

where a = 1 + m* and 0 = m*A ! \ -^ • Note that the critical buckling 

V m* 

loads are independent of the tip body inertia, 

3.4 Stead y -State Solutions of a Cantilevered Beam 
with Tip Body under Axial Thrust (y ^ 0) 

In the previous section, we investigated the buckling of a uni- 
form cantilevered beam with tip body. It will be recalled that in that 
section, the mass center of the tip body was restricted to lie along the 
beam tip tangent line. Presently, we wish to investigate the possi- 
bility of bucklina for the case of nonzero 
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As before, we wish to determine solutions of Eq. (2-19) which are 


indeoendent of time, u = y(x). Ihe governing equation is 


El + P %- [(l - 
,4 ct dx 
dx 


px 
pi. + 


-1 = 

m'^ dx-' 


(3-13) 


dv 

with geometric boundary conditions y(0) = ^ (0) = 0. The natural 
boundary conditions at x = i may be written 


P 

El c cos Y y' ' ' + El y' ’ [1) - me sin y — = 0 (3-14) 

j / pi. + m 


P 

El V"' 01) + ra y' W = 0 

px, + m 


(3-15) 


Integrating Eq. (3-13) and setting the constant of integration to zero, 
we will satisfy the boundary condition Eq. (3-15). 




(1 


px 'I dy 
pi. + m' dx 


0 


(3-16) 


If we make the change of variables 


dx 


w 


and 


z 


P 

r cr /- pi + m -. 2i 1/3|- Px -j 

l-EI p ^ ^ '• pi + m^ 
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the differential equation (3-16) transforms into 


d w 

j 2 

dz 


+ zw = 0 


= 0 goes into z = z^ = f — — 1 ' f- 

1 '■BI ’ '■ 


cr',1/3 rpl + m 2/3 


X = 


X ^ _ m r cr -i1/3 

I goes into z = z^ = ^ —7 — ^ 

2 2/3 ‘•EK p£ + m) ' 

P 


(3-17) 


Substituting Eg. (3-15) into Eq. (3-14) and using the above transforma- 
tion allows us to write the remaining natural boundary condition as 


[(pf + m)Ell^'^^fp pV'^^ 4^ (z ) + me cos y P w(z^) = 

^ cr ' dz 2 cr 2 


-me sin y p 


cr 
(3-18) 


we can write the general solution to Eq. (3-17) as 


" - f=1 ^/3 (| * =2 ■"-1/3 ’i 


'2 3/2 


where c-| and are arbitrary constants. 

Applying the boundary conditions w( z^ ) =0 and Eq. (3-18) we 
obtain the following nonhomogeneous system on c^ and C 2 


" 1 / 3 '“” 


1- J.,/3(«X) 


0 


f/aJ_2^3<3^) + |- /in* c* cos y • X J^^^(PX)lc^ 

* [1 4?- C cos V . X J.,^ 3 ( 6 X> - ^ ^ 2 / 3 lPXll =3 - 

" < 3 -, 9 ) 

where the parameters m*, c* , X, a, P are as defined previously. 
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In general these equations will have a unique solution. If X 
(axial load) assumes a value such that the coefficient matrix in 
Eq. (3-19) IS singular, it can be shown that the equations are inconsis- 
tent. Note that the transcendental equation obtained by setting the 
determinant of the system (3-19) to zero is exactly the same as 
Eq. (3-1 2) except for the replacement c* c* cos y. Hence the system 
(3-19) can become singular for an infinite number of real values of X. 

For those cases where a unique solution exists, we can integrate 
once more and obtain (recall y = 0 at x = 0) 


where 


and 


= _ i rl A- 

£ M a 


■) 


1/3 



3/2 


on the interval 0 < x < £. 

Figure 3-1 illustrates the steady-state solutions given by 
Eq. (3-20) at values of X lying between the first four consecutive cri- 
tical values. Numerical experiments revealed that the overall shapes of 
these curves are insensitive to m*, c*, I, and y within the respective 
ranges of X. 

We have shown that for the case y ^ 0, a unique steady-state 
solution to the vibration equation exists, except when the axial thrust 
assumes certain critical values (such that the linear system (3-19) 
becomes singular) for which no solution exists. 
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CHAPTER 4 


VIBRATION SOLUTIONS 


4 . 1 Natural Frequencies and Mode Shapes (y = 0) — Power 

Series Solution 

In this chapter, we will formulate the eigenvalue problem for the 
transverse vibration of a cantilevered beam with tip body under the ac- 
tion of axial thrust. Since an eigenvalue problem requires the solution 
of a homogeneous differential equation with homogeneous boundary condi- 
tions, we need to restrict y to zero (see boundary condition (2-16). Hie 
a.ial thrust Pq is assumed constant. 

Using Eqs. (2-18) through (2-20), (2-14) and (2-16), we have 


El 


llu 


P I- [(1 - 
o 3x ^ 


px 


pi + 


-) —1 
m"^ 3x ■* 


+ P 


3^u 

3t^ 


= 0 


(4-1) 


u(0,t) 


3u 


0: t > 0 


(4-2) 


I. — , (^,t) + El • c — ^ (£,t) + El • — 4 (£,t) = 0 (4-3) 

3x3t'^ 3x 3x^ 


me 


3 u 
3x3t'^ 


(^,t) - El 


3^ 

3x^ 


(£,t) - m 


pJl + m 3x 


3u „ 3^u 

{i,t) + m ;r 


3t 


(Jl,t) = 0 
(4-4) 
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Assuming a solution 4>(x) , we arrive at the eigenvalue 


problem 


El — + P [(l — 1 - p(D^(j>(x) = 0 (4-5) 

,4 o dx ^ o£ + m ' dx ' ^ 


<j)(0) = (0) = 0 


(4-6) 


2 ^o 


El (j)"'(£) + mfcu + — ■)(!)'(£) + mu) (|>(£) = 0 (4-7) 

^ o£ + m-'^ ^ 


mcP 

2 , , . re- . 2 . o 


mcu(j)(£) + [fl^+mc + £ + ni l ij)' ( - El (|)' ' ( £) = 0 (4-8) 


We wish to find those values of u for which the differential equa- 
tion (4-5), subject to the boundary conditions (4-6) through (4-8) has 
nontrivial solutions. 


Scaling of Eigenvalue problem 

Let 5 = x/£ and define y( ^) s (|)(£5)« Also define the dimension- 
less parameters 


4 p£ 2 

^ El “ 


(4-9) 


m* = m/p£, c* = c/%, j* = 


+ me 
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The eigenvalue problem will then transform into the convenient form 


^ + a^ ^ 




y o) = y' (o) = 0 


y'''(1) + m*(c*b^ + ^ ^ ^ +m*b^y(1) = 


m*c*b'*y(1) + fj^b'* + a^)y* (1) -y"(D = 0 

I + m* 


(4-10) 


(4-11) 


(4-12) 


(4-13) 


Here (') denotes differentiation with respect to 

Let yi(5) and y^ (5) be two eigenfunctions corresponding to 
distinct eigenvalues. It can be shown that the following orthogonality 
condition holds 


1 


/ 


(C)y^ (C)d5 + J*y^(i)y^ (D 


+ m*y^ (l)y^ (1) 


+ m*c*[y^ (l)y^ (1 ) +y^(1)yj(1)] = 0 


cf. Eq. (2-16) of Reference 

The general solution of Eq. (4-10) can be written as a linear com- 
bination of four linearly independent solutions. 

y(?) = ^ ^ ^^4^4^^^ 

For convenience, we specify the initial values of these four functions 
and their first three derivatives as 
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(0) = 6 i = 1,2, 3,4; j = 1,2,3,4 

^ (5. . - Kronecker delta) 

1] 

Applying the geometric boundary conditions (4-11) we find 

c^ = 0 and ^2 ~ ^ 


Application of the natural boundary conditions (4-12) and (4-13) yields 
the homogeneous system of equations on C3 and 04 



where 


3ii = y3"‘(1) + m* (c*b^ + ^ ^ y’(1) + m*b^y2(l) 

2 

^12 " y4’"(1) + m* [c*b + y^d) +m*by^(l) 

(4-14) 

4 / 4 Tn*c* 2 \ 

a^^ = m*c*b y^d) + [j*b + ^ ^ a )y^( 1 ) -y^^'d) 

^22 = n‘*cn\^(1) + (j*b^ a^jy-d) -y^-(l) 


In order for this system to have a nontrivial solution we require 


A(o) ) = 


® 1 1^22 


“ ^ 12^21 


= 0 


(4-15) 


Equation (4-15) determines the natural frequencies (o^. 
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Power Series Expansion for Eigenfunctions 


It still remains to find the functions 73 and 74 . A purely 
numerical technique which is easy to apply consists of integrating two 
initial value problems and iterating on m until values are found which 
satisfy Eq. (4-15). The well-known natural frequencies of a canti- 
levered-free beam serve to bracket the frequencies in the present case of 
a cantilevered beam with a tip body on an accelerating base. We choose 
to address the problem analytically by solving the differential 
Eq. (4-10) in the form of a power series, thus obtaining the eigen- 
functions in functional form. 

Note that for Pq = 0, the differential equation has constant 
coefficients and can therefore be solved in closed form in terms of ele- 
mentary functions (see Reference 1) . For the more general case, Pq * 

0 , the equation is linear with analytic coefficients and has a solution 
which is regular at 5 = 0. The series representing the solution will 
converge for all 5 (see Ince. , E.L., "Ordinary Differential Equations"). 

If we assume a solution of Eq. (4-10) in the form 


insert this expansion and its corresponding derivc*-ives into Eq. (4-10), 
then we obtain the following recursion formulas for the coefficients 


AO O 

= [b Aq + a /(I + m*) A^ - 2a (4-16) 


b'^Aj^ + aVd + m*) • (k + - a^(k + 1 ) (k + 


k+4 


(k + 1) (k + 2) (k + 3) (k + 4) 


(4-17) 
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These conditions determine A., A_, A,, ... once we presribe values for 

4 5 6 

Aq f A^ f A^ f and A^ • 

To generate we must have A^ = 0, A^ = 0, ^2 ~ ^3 

To generate y^(5)f we must have A^ = 0, A^ = 0, A^ = 0, A^ = 1/6 
The eigenfunction Yi(5) corresponding to the eigenvalue is 




4 4 

, o - . 2 ^ a + b? ^ 

, . r1 £.2 a £.4 . a t-5 , 1 t-6 , i 

^12 ^“i^ '■2 ^ ■ 24 ^ 60(1 + m*) ^ 720 ^ *** J 


, . r1 -3 a -3 , 

- ^11<“i’t? ^ -T 20 ^ 


240(1 + m*) 


- a + b. _ 


(4-18) 


E J^ 

AC r then by the ratio test the series 
k=0 ^ 

converges (absolutely) for those values of C for which 


lim 

k+00 



|C| < 1 


We noted above that the differential equation (4-10) has solu- 
tions which are regular at C = 0 and that these series converge for 
all C. Using the recurrence formula ^q. (4-17) we can verify this 
fact. This more elementary procedure was used by Lamb ("Hydro- 
dynamics", p. 335). 

Write the recurrence formula in the form 

\+4 a^ \+2 [a^/(1 + m*) ] (k + 1 ) \+1 

(k + 3) (k + 4) ” (k + 2) (k + 3) (k + 4) A^^ 


(k+1) (k + 2) (k + 3) (k + 4) 
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T&kinq the 


If 


* k +1 ^+2 

lim = L then lim 

k-foo K k-v® ^ 


and 


lim 

k-v® 


^+4 

\ 


limit in the above formula as k>® we obtain the condition 


0 or 


lim 

k-H» 


^+1 


0. Hence the series converges for all E in agreement with 


the general theorv. 


4.2 Free vibration of a Cantilevered Beam with Tip Body 
under Constant Base Acceleration (y ^ 0) 

In this section we investigate the vibrational response of the 
accelerating beam with y ^ 0. It is important to note that in all 
cases considered thus far, y = 0 or steady-state response, a constant 
axial force implied a constant base acceleration, in the present 
case this is no longer true (see Eq. (2-18) . We will assume here that 
a base force Pq( is applied such that the base acceleration a^ 
is constant. 

Neglecting nonlinear terms in u(x,t), the equation of moton is 


KI 



+ a 

o 


a 

ax 


f f p£ + m 



0 


(4-19) 


The boundary conditions are given by Eqs. (2-14), (2-16), and 
(2-20). The boundary condition (2-16) is inhomogeneous. In order to 
solve the differential equation (4-19) we write 


u(x,t) = v(x,t) + f(x) 


(4-20) 


and choose f(x) so that the boundary conditions on v are rendered 
homogeneous. We then have the following requirements on f(x) 


f(0) = 0 


f'(0) = 0 


(4-21 ) 
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0 


(4-21) 


El f " {£) + ma^f ' (£) = 


c El cos Y f ' ' ' (Ji) + El f ' (Ji) = me sin y a 


The partial differential equation on v(x,t) can be written as 


74 ^ ^o ^ + m - px) + p — 

9x 3t 


. {El^t [(pi * „ . px) ^]) 

dx 


(4-22) 


If, in addition to satisfying the conditions (4-21), we require that 
f (x) satisfies 


4 

El ^ + a^ [(pi + m - px) g] =0 (4-23) 

dx 

then the differential equation on v(x,t) , as well as the boundary 
conditions, are rendered homogeneous. 

The conditions on f(x), i.e. , Eq. (4-23) with boundary con- 
ditions (4-21), are precisely those governing the steady-state 
solution of Eq. (4-19) with the boundary conditions (2-14) , (2-16) , 
and (2-20). It was shown in Section 3.4 that f (x) exists for all 
values of ag, except when the system (3-19) becomes singular. f(x) 
IS given by (3-20) with 


X 


3 
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It IS seen that v(x,t) is governed by the system of Equations 


(4-1) through (4-4) with the replacements 


2 2 

c c cos Y, 1^ sin Y 

t t 

P ' p£ + m) 

o r 


thus 


v(x,t) = 2 (a cos 0 ) t + B sin oj t)y ( r) 

- K JC JC iC K 

k=1 

where the frequencies a)]^ are solutions of Eq. (4-15) with 
c* > c* cos y, Po p£ + ™) f the functions Yj^( p) are 

given by Eq. (4-18). 

All frequencies art positive for sufficiently small values of 
3 q. When Sq assimes a critical value such that Eq. (3-12) is sat- 
isfied (with c* -► c* cos ^) , a frequency goes to zero. This cor- 
responds to the situation where f(x) fails to exist. 
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CHAPTER 5 


NATURAL FREQUENCIES AND MODE SHAPES ('f = 0) 
ASSUMED MODES SOLUTION 


In this chapter, approximate solutions to the natural 
frequencies and mode shapes of a cantilevered beam with tip body 
subject to constant axial base acceleration are determined via an 
assumed modes formulation. As indicated in Section 4.1, y must be 
restricted to zero. The mass per unit length, p, and bending 
stiffness, El, are assumed to be constant. 


5.1 Kinetic and Strain Energies 


Let T. , T , and T denote the kinetic energy of the beam, tip 
b t s 

body and total system respectively. Using Eqs. (2-7), (2-8), and 
(2-9) we have in general 





+ (A,t) + c cos Y (^,t)]^} 


T 

s 


1 

2 


(pA + m) 


v^(t) 



me 


’)[■ 


3^u 

'■9x9t 



AM A 

+ me (A,t) [cos T (A,t) - sin Y v^(t)] 

+ j J* (x,t) p dx (5-1) 

o 
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For the case of Y = 0, this reduces to 


= Y (pA + m) v^(t) +1 J* p 


+ {m[|^ (£,t) + c|^ (£,t)]‘ 


(5-2) 


In Appendix A it is shown that the strain energy expression for 
the problem at hand is Eq. (A-11). Using Eq. (2-17) we obtain 


, £ /-2 \2 , £ 
'' ' tt) 

o \3x / o 


(1 + 




(5-3) 


Discretization 


Expressing the beam deformation as 


u(x,t) = ^ (\) {x)q (t) 

i=1 


the energy expro'^-sions (5-2) and (5-3) become 


n n 


i (P^ + m) v^(t) £ ”in^i^i 

° ^ 1=1 ]=1 J 


where 


M 


13 


r 2 

= P f <t>^ (X) <|)^ (X) dx + (I^ + me )(t)^(£)(|)j' (£) + m<i)^ (£)(1)^ (£) 


+ mc[((). (£)(|)^’(£) + <))^(£)(t>^ (£)] 


(5-4) 
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V 


- n n 

= ^ E E 


1=1 j=i 


K. .q.q. 

13 1 3 


where 

I I 

K, = El j* (x)<t.." (X) dx + a^pii j[j - (1 +m*)] (|)’ (x) <)> ’ (x) dx 


- a me (|)! (£)({)'. (X) 


(5-5) 


where 


(') and (•) denote differentiation with respect to x and t, 
respectively 


According to the assumed modes method, the spatial functions 
(x) may be any arbitrary functions which satisfy the geometric 
boundary conditions of the problem and are two times differentiable on 
the interval 0 < x < £. Two sets of functions immediately suggest 
themselves as candidates; the fixed base eigenfunctions of a claraped- 
free beam and those of a clamped beam with tip body. Whereas the 
clamped-free eigenfunctions have the advantage that they can serve for 
beams with any tip body, superior convergence will be realized with 
use of the clamped-tip body eigenfunctions. Furthermore, the latter, 
which were derived in Reference 1, diagonalize the mass matrix, (5-4), 
as a result of their orthogonality property. 

Either set of eigenfunctions satisfy differential equations of 
the form 


= 3%,(x) (5-6) 


4 

For distinct eigenvalues 0^ 
holds 


4 

and 0 j , 


the following indefinite integral 




<t>|^(x)(t)^ (X) 


dx = 


1 


4 4 

(0j - 0^^ 


[0 (b 6' - (!)'>' (j)' ' + 6' ' 4)' ' ' 
i i j 1 j 1 1 
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(5-7) 
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This permits the analytical evaluation of the second integral in Eg. (5-5) 
via integration by parts 


X X 

fij- (1 +m*)] <|)j^ (X) (t)^ (X) dx = f[j- (1 + m*) ] [/(!.’ (x)(j)^(x) dx] dx 


= I — - (1 + ra*) Ife ‘(t.d)! - " (t: ’ + >1)." <t>" ' - efi))' 1 ! 


i (B* - b:> 


1 1 


4 4 

2(6^ + 6p 


£ 4 4 


I ^ - ij>." <1)1 ' + <fc.' <!).''' - 67(1) <l).l 

I 4 _ p4 q>i 9j 9i 9. 


(0^ - B^) y (B. - Bj) 


-2 ef <j). <t). -<!)."' <!)' +<!).’ <i) : "I 

D 1 : 1 1 1 j 


for i ^ 3 . 


Appendix B demonstrates a useful means by which to specialize the 
indefinite integral Eg. (5-7) to the case of i = j 

/[<o:(x)]2dx = f - 24 ).<,v --4-n'‘*’i" 
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BT 

1 


(5-8) 


This permits the evaluation 


f[j - (1 + m*) ] [<{)• ]^ dx = j [| - (1 + m*) ]U[ (<))p 2 - 2<t). 4)- + ^(<1)- ' ) ' 
o ' 3 . 


- <l>" <|)J" + T <t> <t>'- I 

434 " " ^ " Vi 


[(♦;)' - -^7 


' (♦,■•>^*1 ■■ 
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4 '’^i 
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5.2 Eigenvalue Problem 


The application of Lagrange's equations yields the free vibra- 
tion motion equations 


[M]{q} + [K]{q} = { o} 

Seeking harmonic solutions, {q} = {u}e^^^ leads to the eigenvalue 
problem 


([K] - = {0} (5-9) 

The demand of nontrivial solutions of the assumed form yields the set 
of natural frequencies fij- and eigenfunctions 

= ■>'2, n) 

i=1 


The natural frequencies, and eigenfunctions, are the 

sought approximations to those of Section 4. 1 . Note that the values 
of base acceleration for which 

det [K(ao)l = 0 

are approximations to tie critical accelerations of Section 3.3. 
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CHAPTER 6 


FORTRAN PROGRAMS 


Two FORTRAN programs have been created to evaluate the free 
vibration and buckling characteristics of a uniform cantilevered beam 
with tip body under constant axial base acceleration. For each, the 
tip body mass center is restricted to lie along the beam tip tangent 
line. Complete listings of each program, accompanied by annotated 
sample input/output data, are provided in Appendices C and D. 

Hie program of Appendix C computes the natural frequencies and 
mode shapes by implementing the assumed modes formulation of Chapter 
5. Hie eigenfunctions of a clamped-free beam on a fixed base are used 
as the admissible spatial functions. The eigenvalues of the clamped- 

4 

free beam, (S^), are obtained from the roots of the corresponding char- 
ter istic equation 


cos cosh + 1 = 0 

The first fifty of these roots (valid for any beam), appear following 
the sample NAMELIST input data in the appendix. The program functions 
in one of two possible ways, depending upon the values of the param- 
eters AI, AF, and AINC (which denote the initial, final, and incre- 
mental values of the base acceleration, respectively) . If AI = AF 
(AINC arbitrary) , the natural frequencies and mode shapes are computed 
for the single value of base acceleration. If AI < AF, then (only) the 
natural frequencies are computed for the values of base acceleration 
starting with AI and terminating with AF in increments of AINC. Hie 
general algebraic eigenvalue problem Eq. (5-9) is solved via IMSL 
subroutine EIGZS. 
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The program of Appendix D computes the first critical buckling 
load and acceleration by finding the first root of the buckling char- 
acteristic equation. Since available IMSL subroutines evaluate only 
Bessel functions of positive order, the actual characteristic equation 
implemented is not Eq. (3-12) but, rather, an equivalent form expressed 
in terms of positive order Bessel functions 




+ Jl/3(aX) - 2^^ 1 ^) J2/3(BX) + I c*X Js/sCeX) ] = 0 


The IMSL subroutines, MMBSJR and ZBRENT, are used to evaluate the Bes- 
sel functions and the first root of the above characteristic equation, 
respectively. 
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APPENDIX A 


WORK-ENERGY BALANCE 


Starting with the general partial differential equation govern- 
ing transverse beam vibration under axial loading, we derive an 
energy-balance relationship. We then apply the boundary conditions 
ana identify the potential energy expression. This derivation 
obviates the more common geometric arguments which lack a certain de- 
gree of rigor, and reveal the influence of the boundary conditions on 
the strain energy. 

The starting point of our derivation is Eq. (2-4) 



f(x,t) (A-1) 


where we have inserted the term f(x,t) on the right-hand side to ac- 
commodate external loading/length on the beam perpendicular to the x 
axis. The reader can verify that Eq. (A-1) can be written in the al- 
te. nate form 
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3u',2 


3^u>2 
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3x '■3t 3x 


3x 


3x 


3 Fm/ ^ c3u'i2 
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3u 

f(x,t) 


38 


(A- 2) 



This equation holds for 0 < x < il and t > 0. We can then inte- 

0 ^ X ^ 

grate over the region in the xt plane ig < t < t" 


1 r r i'3u'\2 2 /3ui2it . 

2 / [Pllt^ + El(^) 4. t(^) dx 

oX 


r3u 3 3^ui __ 3^u 3^u _ 3u 3u-i£ 

o 3x 3x 


t £ 


t £ 


-i If - //nx,t)|^dxdt 


o o 


o o 


(A-3) 


The right-hand side of Eq. (A-3) represents the work performed by the 
external loading f(x,t) on the beam between t = 0 and t = t. 

The work performed by the axial force Po(t) at the base of 
the beam between t = 0 and t = t is given by 


t 2 t" ^3^ 

r P (t)v (t) dt = - (p£ + m) V (t) 1 - me sin y f (£,t)v dt 

J O O 2 O lt=0 J n 2 O 


■b 3x3t 


(A- 4) 


where we have used Eq. (2-18) (p assumed constant). 

Now the total work done on the system by the external forces f and 
Pq is equal to the change in kinetic plus potential energies. Using 
Eq, (5-1) for the kinetic energy and Eq. (A-3) and (A-4) , this work-energy 
balance can be written as 


I 
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O dX O dXdt 


t .3 


t £ 


- «<'> - i / / f ax .t - 1 (I^ * xc^)[yj ,..X,]^|,^, 


O O 


+ lm [|| l^rt)]^\l^^ + 


[cos Y 1 ^ (£,t) - V^(t) sin y] |^_j 


+ V(t) 



(A- 5) 


V(t) is the potential energy and Q(t) is given by 

tflj(“7r) -- 

Assuming a uniform El and using Eg. (2-12) along with the clamped 
boundary conditions at x = 0 we have 




3x 


2 

8 u T 3u-|£ 

2 3x3t " 3x 3t-lx=0 



Q(^) = / {[ei — I (^ft) + m a^(t) (£,t) 

o Sx"^ 

3 

- me sin Y ■, - (£ft) (£,t)] (£,t) 

Sxat^ 

g2 g2 

- El — I (£,t) (£,t)} dt (A-7) 

9x'^ 
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/ 


The natural boundary conditions (2-14) and (2-16) are eauivalent 


to 


.3 .3 

me cos Y — r (^rt) - El — ^ {l,t) 
3x3t^ 3x-^ 




+ m (Z,t) = 

3t 


(A-8) 


3 2 

(l. + mc^l (£,t) - me cos y a (t) {X.,t) + me cos y (X^,t) 

t ' ^ -.2 O dx ..2 

3 x 3 t 3 t 


+ El — - (£,t) = me sin y • a (t) 

3x^ 


(A-9) 


With the aid of (A-8), (A-9), Eq. (A-7) further simplifies to 
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Q(t) = me cos y “ (Z,t) ^ “ 


3x3t 






- me 


cos y /a^(t) (£,t) (£,t) dt 


3 u 


t 2 

me sin y J ^Q^t) 
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me sin y f 2 (^»t) 

o 3x3t'^ 


dt 
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Now 




3^u 


3x3t 


Usxng the above results, Eq. (A-5) can be put in the form 


V(t) - V(0) =1 / [eI + T(x,t)(|^)2]|^^Q dx 

o 3x 


- me 


t 2 

cos Y /a^(t) -^a,t) dt 


- me 


t a3 
3 u 


sin Y y* 

o 3x3t 


2 / / 3t ^3x^ 
o o 


(A- 10) 


Equation (A-10) gives the potential energy of a clamped beam with tip 
body under axial acceleration agCt) ; El, and p are assumed con- 
stant. Recall that potential energy is defined only up to an arbi- 
trary additive cons*" ant; hence we can set V(0) = 0. 


For the special case y = 0 and constant 
'3x' 


V(t) = { / [ei + T(x) (||)2] dx - Imc a^ ]‘ 


(A-11) 


Admittedly, this derivation of the strain energy does not ap- 
pear to be well motivated and indeed was developed with hindsight. A 
more direct but mathematically sophisticated derivation can be found 
in Reference 2. 
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APPENDIX B 


INDEFINITE INTEGRALS OF EIGENFUNCTIONS 


In this appendix we shall present a method for evaluating the 
indefinite integral 

/[<|)j^(x)J^ dx 

where <|)i(x) is the solution of the differential equation 

.L± = (B-1) 

dx 

corresponding to 6 = Si* This technique is extremely useful in the 
normalization of eigenfunctions and the evaluation of certain inte- 
grals arising in vibration problems in general. One of the earliest 
references to this method is S. Timoshenko ("Vibration Problems in 
Engineering," 2nd edition, p. 335). 

The solution of the differential equation (B-1) depends upon 
the value of 3 as well as x. We indicate this by the notation 
<t> = <J)(x; 3). Formally expanding (|i in a Taylor series we have 

(|)(x; 3 + S3) = <|)(x; 3) + S3 + ... 

Denote by (j»i(x)-the value of <|) corresponding to 3 = and by 
(|)j (x) the value corresponding to 3i + S3. Thus 

(|)^(x) = <|).(x) + (H-)g S3 + ... (B-2) 
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From Eq. (B-1) one can readily show 



4). <!)'.'• - 4). 4)'*' + 4>'4>'' - 4>!4)’' (B-3) 

] 1 1 j 1 j j 1 


where and 0^ are distinct (not necessarily eigenvalues) . 
Substituting Eq. (B-2) into Eq. (B-3) 


- 43^66 


/4>j^(x) [4>^{x) 





f— 1 


5B + 



-4 [4)' • ' 

i '■ i 


dx' 


f— 1 


6B + • . . ] 
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l'- 1 
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where we have used the binomial expansion 


(B-4) 


(B^ + SB)'^ ^ ^i 4B^SB + 


We observe from Eq. (B-1) that the functions 4> depend upon x and B only 
through the argument Bx. Therefore 


34) 3(Bx) _ 34> 

3(Bx) 3B ^3(Bx) 

Similarly 

34> _ 1 34> 

3 (Bx) B 3x 


3 4> 
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Hence 


Ob = '“-=> 

1 1 

Neglecting higher order terms in 66 and using Eq. (B-5) , allows us to wri< e 
Eq. (B-4) as 


M 


dx = <t'.’ <}>V' - ^ 

6 . 1 1 1 6 . ^ 3 1 

1 1 dx 


1 1 H 

+ <ti ' — ■= — (x(!)') - 6'' — — (x(J)') 
^1 6. ^ 2 ' ^i^ ^i 6. dx ' ^i' 
1 dx 1 


Expanding derivatives we obtain the final desired result 




dx = 


X r„4.2 
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. [6 <t> - 2<P’ 4)"’ + (<|)") ] 

4 ■■ 11 1 1 - 1 


467 

1 


— 4) . 4) ! ' ' 
4 ^i ^i 




- 4,; 4," 

4 11 


(B-6) 


Indeed, one can verify this result by direct differentiation. 
Belated integrals, e.g. 



dx 


can be evaluated in the same fashion. 
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APPENDIX C 


NATURAL FREQUENCIES AND MODE SHAPES 
PROGRAM LISTING 


(^ ************************* t ******•****««***•«****«****«*******•********« OQQOQ ^ QQ 

C THIS PROGRAM COMPUTES THE NATURAL FREQUENCIES AND MOOE SHAPES OF A 00000200 
C UNIFORM CANTILEVERED BEAM WITH TIP BODY SUBJECT TO CONSTANT AXIAL 00000300 
C BASE ACCELERATION. AN ASSUMED MODES FORMULATION IS IMPLEMENTED USING 00000400 
C THE EIGENFUNCTIONS OF A CLAMPED-FREE BEAM ON A FIXED BASE AS THE 00000500 
C BASIS THE TIP BODY MASS CENTER IS RESTRICTED TO LIE ALONG THE BEAM 00000600 
C TIP TANGENT LINE. ( WRITTEN BY dOEL STORCH & STEPHEN GATES BASED 00000700 
C UPON CSDL R-1675 OCTOBER 1983 ). 00000800 
C********************************************************************** *00000900 
C 00001000 
(;*********************** input/output files *******»*******»»***********o0001 100 

C 00001200 
C FILE#5. NAMELIST INPUT & "BETA • L" VALUES 00001300 
C FILE#6 OUTPUT ACCELERATIONS S NATURAL FREQUENCIES 00001400 
C FILE/1'7: OUTPUT EIGENFUNCTIONS EVALUATED AT DISCRETE LOCATIONS 00001500 
C 00001600 
C**************** description of namelist input items •**»*»«»***«»*»*»»00001700 


C 



00001800 

c 

"El" 

BEAM BENDING STIFFNESS 

00001900 

c 

"RHO" 

MASS PER UNIT LENGTH OF BEAM 

00002000 

c 


BEAM LENGTH 

00002100 

C 

•MT" 

TIP BODY MASS 

00002200 

C 

•IT" 

TIP BODY INERTIA ABOUT IT'S MASS CENTER 

00002300 

c 

*CT" 

TIP BODY MASS CENTER OFFSET FROM BEAM ATTACHMENT POINT 

00002400 

C 

"AI" 

INITIAL VALUE OF BASE ACCELERATION 

00002500 

C 

"AF" 

FINAL VALUE OF BASE ACCELERATION 

00002600 

c 

•AINC" 

INCREMENT OF BASE ACCELERATION FOR SWEEP 

00002700 

C 

"N* 

NUMBER OF CLAMPED-FREE EIGENFUNCTIONS USED 

00002800 

c 



00002900 


C********************************************************************** *00003000 
C 00003100 
C NOTE ARRAYS DIMENSIONED FOR A MAXIMUM OF 50 DEGREES OF FREEDOM 00003200 
C 00003300 
C********************************************************************** *00003400 




C 


00003500 

ISN 

0002 


IMPLICIT REAL*8 (A-H.O-Z) 

00003600 

ISN 

0003 


REAL*8 L.M.K.MT.IT.MSTAR 

00003700 

ISN 

0004 


DIMENSION K(50,50),M(50,50).BETAL(50).B(50),B4(50).F0L(50). 

00003800 




+F 1L(50), F20(50).F30(50).GAMMA(50. 50). AA( 1275). BB( 1275). D(50), 

00003900 




+Z(50.50).WK(2600).PHI(51.50).Y(51.5) 

00004000 

ISN 

0005 


NAMELIST /INPUT/ El . L . RHO.N.MT. IT .CT . AI . AF . AINC 

00004100 

ISN 

0006 


READ(S.INPUT) 

00004200 

ISN 

0007 


RHOL*RHO*L 

00004400 

ISN 

0008 


MSTAR*MT/RHOL 

00004500 

ISN 

0009 


SRHOL^DSORTCRHOL) 

00004600 



C 


00004700 



C 

REAO BETA»L VALUES 

00004800 



C 


00004900 

ISN 

0010 


DO 1 1*1. N 

00005000 

ISN 

0011 


READ(5.2) BETAL(I) 

00005100 

ISN 

0012 


B(I)*BETAL(I)/L 

00005200 

ISN 

0013 


1 B4(I)=B(I)«*4 

00005300 

ISN 

0014 


2 F0RMAT(E17 8) 

00005400 
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C 



00005500 



C 

ECHO PRINT INPUT DATA 

00005600 



C 



00005700 

ISN 

0015 



WRITE(6,3) 

00005800 

ISN 

0016 


3 

F0RMAT(1H1.24X,'DATA FROM NAMELIST INPUT') 

00005900 

ISN 

0017 



WRITE(6.4) EI.RHO.L 

00006000 

ISN 

0018 


4 

FORMATdHO,' El • ' . 1PE13.6.2X. ' RHO « '.1PE13 6.2X,' L « 

00006100 





+1PE13 6) 

00006200 

ISN 

0019 



WRITE(G.S) MT.IT.CT 

00006209 

ISN 

0020 


5 

FORMATdHO,' MT « '.1PE13 6.2X,' IT * '.1PE13 6.2X,' CT » ' , 

00006218 





+1PE13.6) 

00006227 

ISN 

0021 



WRITE(6,6) AI.AF.AINC 

00006236 

ISN 

0022 


6 

FORMATdHO,' AI « ' . 1PE13.6,2X. ' AF « ' , 1PE13.6,2X, ' AINC » ', 

00006245 





+1PE13.6) 

00006254 

ISN 

0023 



WRITE(6,7) N 

00006263 

ISN 

0024 


7 

FORMATdHO, 'NO. OF CLAMPED-FREE MODES USED • ',12) 

00006272 



C 



00006400 



C 

EVALUATE CLAMPED-FREE EIGENFN & DERIVATIVES P 0 AND L 

00006500 



C 



00006600 

ISN 

0025 



DO 10 I«1,N 

00006700 

ISN 

0026 



BL'BETALd) 

00006800 

ISN 

0027 



S*OSIN(BL) 

00006900 

ISN 

0028 



C*DCOS(BL) 

00007000 

ISN 

0029 



SH=DSINH(BL) 

00007100 

ISN 

0030 



CH*OCOSH(BL) 

00067200 

ISN 

0031 



FOL(D- -2.0D0/SRH0L 

00007300 

ISN 

0032 



F1L(I)=F0L(I)*B(I)«(S*CH + C*SH)/(S»SH) 

00007400 

ISN 

0033 



F20(I)=>F0L(I)*(B(I)**2)*(C+CH)/(S*SH) 

00007500 

ISN 

0034 

10 

F30(I)=F0L(I)*(B(I)«*3)*(S-SH)/(S*SH) 

00007600 



C 



00007610 



c 

ASSEMBLE MASS MATRIX 

00007700 



c 



00007710 

ISN 

0035 



DO 20 1*1, N 

00007800 

ISN 

0036 



DO 20 J>I,N 

00007900 

ISN 

0037 



M(I,d)- (IT+(MT»CT*CT))*F1L(I)*F1L(d) + HT*FOL(I )*FOL(vJ) 

00008000 





++ HT»CT*(F0L(I)«F1L(d) + F1L(I )*FOL(d) ) 

00008 100 

ISN 

0038 

20 

M(J,I)<M(I,J) 

00008200 

ISN 

0039 

- 


DO 21 I«1,N 

00008300 

ISN 

0040 

21 

M(I,I)=M(I,I)+1 ODO 

00008400 



C 



00008410 



C 

ASSEMBLE GAMMA MATRIX 

00008500 



C 



00008510 

ISN 

0041 



DO 30 I«1,N 

00008600 

ISN 

0042 



IP1»I+1 

00008700 

ISN 

0043 



DO 30 d«IP1,N 

00008800 

ISN 

0044 



BMB>B4(I)-B4(0) 

00008900 

ISN 

0045 



T1-(1 000+MSTAR)*(F20(I)*F30(d) - F30( I )*F20(d) ) 

00009000 

ISN 

0046 



T2=MSTAR*(B4(I)*F0L(I)‘F1L(d) - B4(d)»F1L(I)»F0L(d)) 

00009100 

ISN 

0047 



T3»2 0D0«(B4(I)+B4(d))*F20(I)*F20(d)/(L*BMB) 

00009200 

ISN 

0048 



T4»4.O0O*B4(I)«B4(U)*F0L(I)»FOL(d)/(L*BMB) 

00009300 

ISN 

0049 



GAMMA(I,d)»(T1-T2-T3+T4)/BMB 

00009400 

ISN 

0050 

30 

GAMMA(J,I)«GAMMA(I,U) 

00009500 

ISN 

0051 



DO 40 1*1, N 

00009600 

ISN 

0052 



T1«MSTAR*(L*(F1L(I)**2) ♦ 3 000»FOL(I )*F1L( I ) )/4 ODO 

00009700 

ISN 

0053 



T2*d 0DO+MSTAR)*F20(I)*F3O(I)/(4 0D0*B4(I)) 

00009800 

ISN 

0054 



T3*(L*F1L(I))»*2 + 3 ODO*(FOL( I )»*2 ) ♦ ( F20( I ) ••2 )/B4( I ) 

00009900 

ISN 

0055 

40 

GAMHA(I,I)--T1-T2-T3/(8 ODO»L) 

00010000 



C 



00010010 
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C ASSEMBLE STIFFNESS MATRIX 

00010100 



C 


00010110 

ISN 

0056 


AX^AI 

00010200 

ISN 

0057 

50 

DO 60 I>1,N 

00010300 

ISN 

0058 


IP1«I+1 

00010400 

ISN 

0059 


00 60 d»IP1.N 

00010500 

ISN 

0060 


K(I,d)«AX RHaL«GAMMA(I.d)-MT*CT»AX‘F1L(I)*F1L(d) 

00010600 

ISN 

0061 

60 

K(J.I)-K( ,J) 

00010700 

ISN 

0062 


DO 70 I»1 N 

00010800 

ISN 

0063 

70 

K(I.I)=EI*B4(I)/RH0 + AX*RH0L*GAMMA(I.I)-HT»CT*AX»(F1L(I)**2) 

00010900 



C 


00010910 



C LOAD K INTO AA & H INTO BB BOTH SYMMETRIC STORAGE 

00011000 



C 


00011010 

ISN 

0064 


II»1 

00011100 

ISN 

0065 


DO 80 1*1, N 

00011200 

ISN 

0066 


DO 80 0*1,1 

00011300 

ISN 

0067 


AA(II)»K(I,d) 

00011400 

ISN 

0068 


BB(II)*H(I,d) 

00011500 

ISN 

0069 


11=11+1 

00011600 

ISN 

0070 

80 

CONTINUE 

00011700 



C 


00011710 



C COMPUTE EIGENVALUES S EIGENVECTORS 

00011800 



C 


00011810 

ISN 

0071 


IdOB=2 

00011900 

ISN 

0072 


IZ»50 

00012000 

ISN 

0073 


CALL EIGZSCAA.BB.rr.IOOB.D.Z.IZ.WK.IER) 

00012100 



C 


00012101 



C OUTPUT FIRST 5 EIGENVALUES & ACCELERATION 

00012200 



C 


00012210 

ISN 

0074 


WRITE(6,90 AX 

00012300 

ISN 

0075 

90 

FORMAT! 1H0, 'BASE ACCELERATION • '.1PE13.6) 

00012400 

ISN 

0076 


WRITE(6,91 ) 

00012500 

ISN 

0077 

91 

FORMAT! 1H0 4X, 'NATURAL FREQUENCIES') 

00012600 

ISN 

0078 


WRITE!6,92, 

00012610 

ISN 

0079 

92 

FORMAT! 1H ,1X.'N'.5X.'0MEGA!N)**2 !SEC**-2)') 

00012620 

ISN 

0080 

- 

DO 93 1-1,5 

00012630 

ISN 

0081 

93 

WRITE!6,94) I,D!I) 

00012640 

ISN 

0082 

94 

F0RMAT!1H I2,5X,1PE13 6) 

00012650 



C 


00012700 



C TEST FOR EIGENVECTOR COMPUTATION OR ACCEL SWEEP 

00012800 



C 


00012810 

ISN 

0083 


IF! AI .EO AF ) GO TO 100 

00012900 

ISN 

0085 


IF! AX .GE. AF ) GO TO 1000 

00013000 

ISN 

0087 


AX»AX+AINC 

00013100 

ISN 

0088 


GO TO 50 

00013200 



C 


00013210 



C EVALUATE CLAMPL3-FREE EIGENFUNCTIONS AT 51 STATIONS 

00013300 



C 


00013310 

ISN 

0089 

100 

0X*L/5O 000 

00013400 

ISN 

0090 


00 110 1*1, N 

00013500 

ISN 

0091 


BL»BETAL!I) 

00013600 

ISN 

0092 


SS=0SIN!BL)«0SINH!BL) 

00013700 

ISN 

0093 


SMS«OSIN!BL) - 0SINH!BL) 

00013800 

ISN 

0094 


CPC*DC0S!BL) + 0C0SH!BL) 

00013900 

ISN 

0095 


X=O.ODO 

00014000 

ISN 

0096 


DO 105 0*1,51 

00014100 

ISN 

0097 


bx=b!i)*x 

00014200 

ISN 

0098 


phi!o,i)-sms*!osin!bx)-dsinh!bx)) + cpc*!dcos!bx)-dcosh!bx)) 

00014300 


48 


ISN 

0099 


PHI(d.I)-PHI(d,I)/(SRHOL*SS) 

00014400 

ISN 

0100 

105 

X»X+OX 

00014500 

ISN 

0101 

110 

CONTINUE 

00014600 



C 


00014610 



C COMPUTE FIRST 5 EIGENVECTORS 

00014700 



C 


00014710 

ISN 

0102 


00 125 IP«1.51 

00014800 

ISN 

0103 


00 120 IM»1,5 

00014900 

ISN 

0104 


SUM=0 000 

00015000 

ISN 

0105 


DO 115 d»1,N 

00015100 

ISN 

0106 


PROD»PHI(IP.d)*Z(d,IM) 

00015200 

ISN 

0107 

115 

SUH*SUM+PROO 

00015300' 

ISN 

0108 

120 

Y(IP,IH)*SUH 

00015400 

ISN 

0109 

125 

CONTINUE 

00015500 



C 


00015510 



C OUTPUT EIGENVECTORS FOR PLOTTING 

00015600 



C 


00015610 

ISN 

0110 


X»0 ODO 

00015700 

ISN 

0111 


DO 130 1-1,51 

00015800 

ISN 

0112 


WRITE(7,131) X.(Y(I,d),d-1.5) 

00015900 

ISN 

0113 

130 

X-X+DX 

00016000 

ISN 

0114 

131 

F0RMAT(6(1X,E12 5)) 

00016100 

ISN 

0115 

1000 

STOP 

00016200 

ISN 

0116 


END 

00016300 
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EXAMPLE PROBLEM PARAMETERS 


BEAM PARAMETERS 
7 2 

El = 6 25 • 10 Ib-ft 
L = 1340 ft 
p = 0 41 72 slug/ft 


TIP BODY PARAMETERS 
m = 1 55 28 slugs 
Ij = 0 0 slug-ft^ 

. » 00ft 


NAMELIST INPUT AND CLAMPED-FREE BEAM "EIGENVALUES" 


&INPUT EI=6 25E7,L»134 O.RHO* .4172 ,N=20 
MT«1S5.28,IT*0 0,CT=0.0 
AI=0 0,AF»60 0.AINC*20 O.&END 


18t51041E+01 
4694091 1E+01 
78547574E+01 
10995541E+02 
141371G8E+02 
17278760E+02 
20420352E+02 
23561945E+02 
0.26703538E+02 
0 29845 130E+02 
0 32986723E+02 
0 36128316E+02 
0.392G9908E+02 
0.42411501E+02 
0 45553093E+02 
0 48694686E+02 
0 51836279E+02 
0.54977871E+02 
0 58119464E+02 
61261057E+02 
64402649E+02 
67544242E+02 
70685835E+02 
73827427E+02 
7G969020E+02 
80110613E+02 
83252205E+02 
86393798E+02 
8953539 1E+02 
92G76983E+02 
95818576E+02 
989G01G9E-»02 
1021017GE+03 
10524335E+03 
10838495E+03 
11152G54E+03 
11466813E+03 
11780972E+03 
12095132E+03 
12409291E+03 
12723450E+03 
0. 130376 10E+03 
0. 133517G9E+03 
13665928E+03 
13980087E+03 
14294247E+03 
14608406E+03 
14922565E+03 
15236724E+03 
15550884E+03 


\ 


FIRST FIFTY ROOTS iffi) OF 

cos cosh (3 2 = -1 
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PROGRAM OUTPUT 
ACCELERATION SWEEP RUN 




DATA 

FROM NAMELIST 

INPUT 



El • 

S.250000E+07 

RHO « 

4 172000E-01 

L » 

1 

340000E+02 

MT » 

1.552800E+02 

IT • 

0 OOOOOOE+00 

CT • 

0 

OOOOOOE+00 

AI = 

0 OOOOOOE+00 

AF • 

6 OOOOOOE+01 

AINC • 

2 

OOOOOOE+01 


NO OF CLAMPED-FREE MOOES USED • 20 


BASE ACCELERATION • 0 OOOOOOE+00 

NATURAL FREQUENCIES 
N 0MEGA(N)»*2 (SEC**-2) 

1 4.625102E-01 

2 1 151594E+02 

3 1 176526E+03 

4 5 085561E+03 

5 1 483397E+04 

BASE ACCELERATION • 2.000000E+01 

NATURAL FREQUENCIES 
N 0MEGA(N)*»2 (SEC**-2) 

1 2 785320E-01 

2 1 09411GE+02 

3 1 155G8GE+03 

4 5 039933E+03 

5 1 475390E+04 

BASE ACCELERATION » 4 OOOOOOE+01 

NATURAL FREQUENCIES 
N 0MEGA(N)»*2 (SEC»*-2) 

1 9.289338E-02 

2 1 03G581E+02 

3 1 134843E+03 

4 4 994303E+03 

5 1 4G7382E+04 

BASE ACCELERATION « G OOOOOOE+01 

NATURAL FREQUENCIES 
N 0MEGA(N)»»2 (SEC«*-2) 

1 -9 4GG231E-02 

2 9.789877E+01 

3 1 113999E+03 

4 4 948672E+03 

5 1 459374E+04 
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PROGRAM OUTPUT 
MODE SHAPES RUN 


DATA FROM NAMELIST INPUT 


El • 

6 

250000E+07 

RHO • 

4 

172000E-01 

L • 

1 

340000E+02 

MT ■ 

1 

552800E+02 

IT - 

0 

OOOOOOE+00 

CT - 

0 

OOOOOOE+00 

AI « 

3 

220000E+01 

AF ■ 

3 

220000E+OI 

AINC • 

0 

OOOOOOE+00 


NO OF CLAMPEO-FREE MOOES USED • 10 

BASE ACCELERATION • 3 220000E«0I 

NATURAL FREQUENCIES 
N 0MEGA(N)**2 (SEC**-2) 

« 1 6S5051E-0t 

2 1 0S9027E*02 

3 1 I429T2E*03 

4 5 0I2099E«03 

5 1 47050SE*04 


EIGENFUNCTIONS (NON-NORMAUZED) 


X 

V(x) 

v<2) 

''(x) 

y(3) 

y(4) 

''(x) 

v<5) 

''(x) 

0 00000 E* 00 

0 OOOOOE+00 

0 OOOOOE+00 

0 OOOOOE+00 

0 OOOOOE+00 

0 OOOOOE+00 

0 26800E«01 

0.41721E-04 

-0 80237E-03 

0 25418E>02 

•0 52092E-02 

0 86654E-02 

0 53600E*01 

0. 16641E-03 

•0 31241E-02 

0 96775E-02 

-0 19296E-01 

0 31340E-01 

0 60400E401 

0 37331E-03 

•0 68381E’02 

0 20664E‘01 

*0 39984E-01 

0 83036E*01 

0 10720e*02 

0 66154C-03 

•0 f18t9E-Ot 

0 34748E-01 

•0 65064E-01 

0 9BB44E-01 

0 13400C+02 

0 10302E-02 

•0 17939E-01 

0 51180E-01 

•0 92400E-01 

0 13415E*00 

0 16080e*02 

0 14782E-02 

-0.25072E-01 

0 69230E-01 

-0.11993E+00 

0 16493E+00 

0 18760E^02 

0 20047E-02 

-0 33086E-0I 

0 88199E-01 

“0 14574E+00 

0 18787E+00 

0 21440E+02 

0 2608SE-02 

-0 4i853E-OI 

0 I074U400 

-0. iG814E+00 

0 2004SE+00 

0 24120E+02 

0 3288SE-02 

-0 51248E-01 

0 12620E+00 

•0.18570E*00 

0 20108E+00 

0 26800E+02 

0 4043SE-02 

-0 6U45E-01 

0 t4395E+00 

•0 19733E+00 

0 109O8E+OO 

0 29480E+02 

0 48722E-02 

-0 714t9E-Of 

0 I6012E*00 

•0 20219E+00 

0 1S494E+00 

0 32160E+02 

0 57733E-02 

•0 8t947E-OI 

0 I7422E*00 

-0 19980E+00 

0 13012E*00 

0 34840E+02 

0 67456E-02 

•0 92607E-01 

0 18581E400 

•0 19001E+00 

0 86888E'01 

0 37520E*02 

0 77875E-02 

•0. 10328E*00 

0 19452E*00 

-0 17307E*00 

0 38085E-01 

O 40200C*02 

O 88976£-02 

•O 11386C*00 

O 20006t*00 

-0 149596*00 

-O 130196-01 

0 42880E+02 

0 10074E-01 

•0 12424E*00 

0 20222E+00 

-0 12041E+00 

-0 G2B8GE-01 

0 45560E+02 

0 11317E-0t 

■0 13430E*00 

0 2OO80E+OO 

•0 86S20E-01 

-0 10798E*00 

0 48240E*02 

0 12622E-01 

•0 14394E*00 

0 19600E400 

•0 49540E-01 

-0 14514E*00 

0 50920E+02 

0 13990E-01 

■0 15308E*00 

0 18790E400 

•0 10G76E-01 

-0 17180E*00 

0 53600E+02 

0 15418E-01 

0 18ISIE*0O 

0 17579E400 

0 28371E-01 

-0 18620E+00 

0 56280E*02 

0 16905E-01 

•0 16947E*00 

0 16076E400 

0 66022E-01 

-0 18730E*00 

0 58960E+02 

0 18449E-01 

•0 178S6E*00 

0 14279E*00 

0 10078E*00 

-0 17499E*00 

0 61940E«02 

0 20048E-01 

-0 18283E*00 

0 *2219E400 

0 13123E*00 

-0 15012E*00 

0 64320E+02 

0 21700E-01 

-0 18819E*00 

0 9VJO0E-01 

0 i5610E+00 

-0 11453E+00 

0 67000E*02 

0 23404E-01 

■0 1929!E*00 

0 74706E-01 

0 17434E*00 

-0 70829E-01 

0 69680E^02 

0 25158E-01 

■0 19602E*00 

0 48708E-01 

0 18521E*00 

-0 22121E-01 

0 72360E+02 

0 26960E-01 

•0 19840E*00 

0 21882E-01 

0 18832E*00 

0 282tSE-01 

0 75040E*02 

0 280O8E-O1 

•0 19989E*00 

•0 524tBE-02 

0 18356E^00 

0 76647E-01 

0 77720E+02 

0 30700E-01 

•0 19989E+00 

•0 32I43E-0I 

0 171ltE*00 

0 11965E+00 

0 80400E+02 

0 32634E-01 

•0 19897E+00 

-0 58310E'01 

0 15I44E*00 

0 15405E*00 

0 83080E+02 

0 34608E-01 

-0 19693E*00 

•0 83221E-01 

0 12541E*00 

0 17745E+00 

0 85760E*02 

0 36620E-01 

-0 19375E*00 

-0 10634E400 

0 94224E-01 

0 1884€E*00 

0 88440E+02 

0 38€69E’01 

-0 18945E*00 

-0 12725E+00 

0 59078E-01 

0 18805E*00 

0 91 120E+02 

0 40752E-01 

-0 18405E+00 

•0 14555E*00 

0 21375E-01 

0 17023E+00 

0 93800E+02 

0 42888E’01 

-0 1775SE+00 

-0 18086E+00 

-0 17201E-01 

0 14236E*00 

0 96480£*02 

0 45014E-0t 

•0 17001E+00 

-0 17282E4QO 

-0 54938E-01 

0 t0464E*00 

0 99160E*02 

0 47188E-01 

•0 16I44E*00 

•0 I8125E+00 

■0 90385E-01 

0 59478E-Ot 

0 10184E+03 

0 49388E'01 

•0 15190E*00 

•0 I8604E+00 

-0 12225E+00 

0 96639E-02 

0 10452E*03 

0 S16I2E<01 

•0 14143E+00 

•0 I8701E*00 

*0 14896E+00 

-0 40771E-01 

0 10720E*03 

0 53859E-01 

•0 13009E»00 

•0 184I0E400 

-0 16925E*00 

-0 87855E-01 

0 10988E+03 

0 56126E-01 

-0 11792E*00 

•0 17739E+00 

'0 18237E+00 

-0 12853E»00 

0 11256E+03 

0 58411E-01 

-0 10502E+00 

-0 16705E400 

-0 18787E+00 

-0 16008E*00 

0 11524E*03 

0 60712E-01 

-0 91462E-0! 

•0 I5338E^00 

-0 18579E*00 

-0 18089E*00 

0 11792E+03 

0 63026E-0! 

-0 77315E-01 

-0 13660E*00 

-0 17603E*00 

•0 189ltE+00 

0 12060E+03 

0 65353E'01 

-0 62637E-01 

•0 I1890E+00 

■0 15865E*00 

-0 18335E*00 

0 12328E*03 

0 67690E-01 

-0 47501E'01 

-0 94618E-01 

•0 13414E*00 

-0 16359E+00 

0 12596E403 

0 70036E-01 

-0 32000E-01 

-0 70237E-01 

•0 10368E+00 

-0 13153E+00 

0 12864E«03 

0 72387E-01 

-0 16262E-01 

•0 44470E-01 

•0 69197E-01 

-0 90955E-0I 

0 13132E+03 

0 74741E-01 

-0 37928E-03 

-0 17859E-01 

•0 32217E-01 

•0 44939E-01 

0 13400E+03 

0 77096E-01 

0 15543E-01 

c 

0 899I4E-02 
zn 

0 54750E-02 

0 26399E-02 
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APPENDIX D 


«* 


BUCKLING LOAD/ACCELERATION 
PROGRAM LISTING 


(^ «**•««•«•*«•«***«•*»«********«***•*****••*»******•********* t *«*•****«*) QQQQQ ^ QQ 

C THIS PROGRAM COMPUTES THE FIRST CRITICAL BUCKLING LOAD/ACCELERATION 00000200 
C OF A UNIFORM CANTILEVERED BEAM WITH TIP BODY THE "BUCKLING LOAD" 00000300 


C CORRESPONDS TO THE BEAM ROOT AXIAL FORCE AND IS DIRECTLY PROPOR- 00000400 
C TIONAL TO THE "BUCKLING ACCELERATION". THE TIP BODY HASS CENTER IS 00000409 
C RESTRICTED TO LIE ALONG THE BEAM TIP TANGENT LINE ( WRITTEN BY 00000418 
C JOEL STORCH S STEPHEN GATES BASED UPON CSDL R-1675 OCTOBER 1983 ) 00000427 

(;*«******«*«*«****•****•*«*******•«***«******«*«•*** •***************»***QQQQg^30 

C 00000445 
Q********************** INPUT/OUTPUT FILES ••»»»»»»»»»»»»»»»»»**»»»»»»*»00000454 
C 000004G3 
C FILE#5: INPUT NAMELIST INPUT DATA 00000472 
C FILE#6: OUTPUT CRITICAL BUCKLING LOAD/ACCELERATION 00000481 
C 00000490 
C«******«******* description OF NAMELIST INPUT ITEMS ••***»»»»*»«***»*»»00000500 
C' 00000509 
C "El" BEAM BENDING STIFFNESS 00000518 
C "RHO" MASS PER UNIT LENGTH OF BEAM 00000527 
C "L" BEAM LENGTH 00000536 
C "MT* MASS OF TIP BODY 00000545 
C "CT" TIP BODY MASS CENTER OFFSET FROM BEAM ATTACHMENT POINT 00000554 
C 00000563 




(^*«* ««**«««***«*«««*«*«*«««««♦«««««*««««**«**«******♦••♦****** **********f()0QQ572 



c 



C 000 1000 

ISN 

0002 


IMPLICIT REAL*8(A-H.0-Z) 


00001100 

ISN 

0003 


REAL«8 L.LAMS.MSTAR.MT 


C0001200 

ISN 

0004 


DIMENSION LAMS(2) 


00001300 

ISN 

0005 


DATA LAMS/1 00-8.1 87D0/ 


00001400 

ISN 

0006 


COMMON ALPHA, BETA. Cl, C2.C3. 01. 02 


00001500 

ISN 

0007 


EXTERNAL F 


00001600 

ISN 

0008 


NAMELIST /INPUT/ El . L.RHO.MT.CT 


(0001610 

ISN 

0009 


REA0(5, INPUT) 


00001620 

ISN 

0010 


MSTAR*MT/(RHO»L) 


00001630 

ISN 

0011 


CSTAR»CT/L 


00001640 



c 



00001650 



c 

ECHO PRINT INPUT DATA 


00001653 



c 



00001656 

ISN 

0012 


WRITE(6,1) El, RHO, L 


00001660 

ISN 

0013 


1 F0RMAT(1H1,24X.'0ATA FROM NAMELIST INPUT' 

.//.' El » 

',1PE13.6,2X 00001670 




+,' RHO • '.1PE13.6.2X,' L • '.1PE13.6) 


C >001680 

ISN 

0014 


WRITE(6,2) MT.CT 


0C0O1690 

ISN 

0015 


2 FORMAT(1HO,' MT • '.1PE13 6.2X.' CT » ' 

, 1PE13.6, ) 

00001700 



c 



00001710 



c 

COMPUTE CONSTANTS 


00001720 



c 



00001730 

ISN 

0016 


ALPHA* 1 000+MSTAR 


00002000 

ISN 

0017 


BETA»MSTAR/0S0RT(1 000+1 000/MSTAR) 


00002100 

ISN 

0018 


C1*4 000/(3 ODO"ALPHA) 


00002200 

ISN 

0019 


C2*0S0RT( ALPHA) 


00002300 

ISN 

0020 


C3*2.000»C2/(3 000*BETA) 


00002400 

ISN 

0021 


C4«1 5DO*DSORT(MSTAR) 


00002500 

ISN 

0022 


C5*2 O0O*0S0RT(MSTAR)/BETA 


00002600 

ISN 

0023 


01=C4*CSTAR 


00002900 
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ISN 

0024 



D2=C2-C5»CSTAR 


00003000 

ISN 

0025 



A’LAMS(I) 


00003200 

ISN 

0026 



B«LAMS(2) 


00003300 

ISN 

0027 



MAXFN«80 


00003400 



C 




00003410 



C 

COMPUTE FIRST ROOT OF BUCKLING CHARACTERISTIC 

EQUATION 

00003420 



C 




00003430 

ISN 

0028 



CALL ZBRENT(F,1 00-8,S.A,B,MAXFN,IER) 


00003500 

ISN 

0029 



IF(IER EO. 0) GO TO 4 


00003600 

ISN 

0031 



WRITE(6,3) lER 


00003700 

ISN 

0032 


3 

FORMAT(1HO,')RROR IN ZBRENT ALGORITHM lER 

• ',13) 

00003800 

ISN 

0033 



STOP 


00004000 

ISN 

0034 


4 

PCR = (2.25DO*f /L“2)»B**2 


00004100 

ISN 

0035 



ACR=PCR/(RHO» + MT) 


00004200 

ISN 

0036 



WRITE(6,5)PCR,ACR 


00004300 

ISN 

0037 


5 

FORMATdHO, 'FIRST CRITICAL BUCKLING LOAD » 

' , 1PE13.6,/, 

00004400 





♦'FIRST CRITICAL BUCKLING ACCELERATION • ', 

1PE13 6) 

00004500 

ISN 

0038 



STOP 


00005400 

ISN 

0039 

- 


END 


00005500 




Qtt*^************** **♦**♦****♦*♦♦•****•**•***♦♦♦•*« ***♦♦*♦*•♦♦*♦♦****♦♦000055^0 



C 

THIS FUNCTION SUBPROGRAM EVALUATES THE BUCKLING CHARACTERISTIC 

00005520 



C 

EQUATION. EO (3-12) EXPRESSED IN TERMS DF BESSEL FUNCTIONS OF 

00005530 



C 

POSITIVE ORDER. 

00005540 



(^««««««****««*««*«*«««««*««*««*««**«***«****** *******«*♦« **««*««*«****«*QQQQ555Q 

ISN 

0002 ‘ 


FUNCTION F(LAM) 

00005600 

ISN 

0003 


IMPLICIT REAL*8(A-H,0-Z) 

00005700 

ISN 

0004 


REAL*8 LAM 

00005800 

ISN 

0005 


COMMON ALPHA, BETA, Cl. C2.C3.D1.D2 

00005900 

ISN 

0006 


DIMENSION T(8),R0(2).WK(4) 

00006000 

ISN 

0007 


A1*ALPHA»LAM 

00006100 

ISN 

0008 


A2*BETA*LAM 

00006200 



C 


00006300 



c 

COMPUTE BESSEL F, 'NOTIONS OF ORDERS 1/3. 2/3, 4/3, 5/3 

00006400 



c 

AT ARGUEMENTS A1 AND A2 

00006500 



c 


00006600 

ISN 

0009 


IS«0 

00006700 

ISN 

0010 


DO 40 1=1.2 

00006800 

ISN 

0011 


OR*OFLOAT(I)/3.0DO 

00006900 

ISN 

0012 


DO 30 NA=1,2 

00007000 

ISN 

0013 


IF(NA EO 1) ARG=A1 

00007100 

ISN 

0015 


IF(NA EO 2) ARG=A2 

00007200 

ISN 

0017 


CALL MM8SJR(ARG,0R,2.RU.WK.IER) 

00007300 

ISN 

0018 


IF( lER EO 0) GO TO 20 

00007400 

ISN 

0020 


WRITE(6,10) lER.ARG.OR 

00007500 

ISN 

0021 

10 

FORMATdHO, 5X, 'ERROR IN EVALUATING BESSEL FUNCTION' . 3X ,' IER = ' . 

00007600 




+ I3,2X,'ARGUEMLVT*'.E13.5,2X.'0RDER»'.F7 5) 

00007700 

ISN 

0022 


STOP 

00007800 

ISN 

0023 

20 

IS»IS*1 

00007900 

ISN 

0024 


T(IS)»RJ(1) 

00008000 

ISN 

0025 


IS=IS+1 

00008100 

ISN 

0026 


T(IS)»RvJ(2) 

00008200 

ISN 

0027 

30 

CONTINUE 

00008300 

ISN 

0028 

40 

CONTINUE 

00008400 



c 


00008500 



C 1 

COMPUTE FUNCTION "F" 

00008600 



C 


00008700 

ISN 

0029 


F=(C1/LAM*T(5)-T(6))*(C3/LAM*T(3)-C2*T(4)+D1«LAM*T(3))+ 

00008800 




1 T(1)*(02*T(7)+01*LAM*T(8)) 

00008900 

ISN 

0030 


RETURN 

00009000 

ISN 

0031 


END 

00009100 
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SAMPLE INPUT/OUTPUT FOR BUCKLING PROGRAM 


BEAM PARAMETERS 
El =■ 6 25 • 10’lb-ft^ 
L - 134.0 ft 
p =■ 0 41 72 slug/ft 


TIP BODY PARAMETERS 

m = 155 23 slugs 
c - 00ft 


NAMELIST INPUT 

&INPUT EI»S.25E7,RH0« 4172.L«134 0.MT«155.28,CT*00.0,&EN0 


PROGRAM OUTPUT 
DATA FROM NAMELIST INPUT 

El • 6 250000E+07 RHO • 4 172000E-01 L » 1 340000E+02 

MT • 1 552800E+02 CT • O.OOOOOOE+00 
FIRST CRITICAL BUCKLING LOAD » 1 05450SE+04 lb 

2 

FIRST CRITICAL BUCKLING ACCELERATION • 4 993284E+01 ft/s 

BEAM PARAMETERS TIP BODY PARAMETERS 

SAME AS ABOVE m = 155 28 slugs 

c =■ 67 0 ft 


PROGRAM OUTPUT 

DATA FROM NAMELIST INPUT 

El • 6 250000E+07 RHO » 4 172000E-01 L • 1.340000E+02 

MT » 1 552800E+02 CT » 6.70OC0OE+O1 

FIRST CRITICAL BUCKLING LOAD • 5 261190E+03 lb 

FIRST C''ITICAL BUCKLING ACCELERATION • 2.491273E+01 ft/s^ 
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APPENDIX E 


COMPARISON OF SELECTED NUMERICAL Rj SULTS 


In this appendix we compare the numerical values for the funda- 
mental natural frequencies obtained from four distinct numerical 
procedures. Uie results are for the case of a tij body with mass 
center alonq the beam tip tangent line and the particular parameters 
listed below. 

Beam Parameters Tip Body Parameters 

El = 6.25 • 1q 7 slug-ft2 m = 155. 2J slug 

=» 134.0 ft It = 5.0 1q 5 slug-ft2 

p = 0.4172 slug/ft c = 67.0 rt 

The first two columns of natural frequencies are based upon the 
assumed modes formulation of Chapter 5, using the first 10 fixed base 
eigenfunctions of a clamped-free beam and clamped beam with tip body, 
respectively. The spline based .Galerkin method used to generate the 
results of the third column, is described in Referc^ce 2. The results 
of the fourth column are based upon the power series solution of Sec- 
tion 4.1, using six terms in each series. 
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t An t 


Fundamental Natural Frequencies (s' 

-2) 

Axial Base 
Acceleration 
(ft/s2) 

Assumed-Modes Method 
10 Clamped — Free 
Eigenfunctions 

Assumed-Modes Method 
10 Clamped — Tip Body 
Eigenfunctions 

Spline-Based 
Galerkin Method 
1 0-cubic Splines 

Power-Series 

Solution 

6-Terms in Each Series 

0.0 

0.1329 

0.1312 

0.1312 

0.1312 

2.4913 

0.1198 

0.1181 

0.1181 

0.1181 

4.9825 

0.1067 

0.1050 

0.1050 

0.1050 

7.4738 

0.0936 

0.0919 

0.0919 

0.0919 

9.9651 

0.0805 

0.0788 

0.0788 

0.0788 

12.4564 

0.0674 

0.0657 

0.0657 

0.0657 

14.9476 

0.0543 

0.0526 

0.0526 

0.0526 

17.4389 

0.0411 

0.0395 

0.0395 

0.0395 

19.9302 

0.0280 

0.0263 

0.0263 

0.0263 

22.4215 

0.0148 

0.0132 

0.0132 

0.0132 

24.9127 

0.0017 

0 

0 

0 
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